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Abstract. A modified version of the Faddeev three-body equation to accom- 
modate the Coulomb interaction, which was used in the study of three-nucleon 
^sO | bound states, is applied to the proton-deuteron scattering problem at energies 

below the three-body breakup threshold. A formal derivation of the equation 
■ in a time-independent scattering theory is given. Numerical results for phase 

shift parameters are presented to be compared with those of another methods 
and results of the phase shift analysis. Differential cross sections and nucleon 
analyzing powers are calculated with the effects of three-nucleon forces, and 
^ | these results are compared with recent experimental data. The difference be- 

tween the nucleon analyzing power in proton-deuteron scattering and that in 
neutron-deuteron scattering is discussed. 

X 

1 Introduction 

Recently, precise measurements of observables in proton-deuteron (pd) elastic 
scattering have been performed at several energies [||, ^ ||, ||, H), ^, 0, §| ||, [u], |llj 
for the details of nuclear forces, such as an off-shell difference in realistic nucleon- 
nucleon (NN) potentials or the evidence for three-nucleon forces (3NFs). For 
obtaining less ambiguous information from these experimental data, we need 
to calculate the corresponding observables accurately by solving the problem of 
the three-body system consisting of two protons and one neutron. As is well 
known, however, it is not trivial task to treat this three-body system because of 



the long-range Coulomb force acting between two protons. In ref. [12|, Sasakawa 
and Sawada proposed a modification of the three-body Faddeev equation |^] to 
accommodate Coulomb forces. The Coulomb- modified Faddeev equations were 



solved in refs. fl4| , 15, T^] to obtain the wave function of 3 He bound state for real- 
istic NN potentials with or without 3NFs, and were shown to work well both in its 
accuracy and the speed of the computation. In this paper, we apply the method 
to the pd scattering problem and present some numerical results. The formulary 
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aspect of the method in refs. [15, 16] is that we express the Faddeev equation as 



an integral equation in coordinate space, and then, solve the equation with an it- 



erative method, which we called the Method of Continued Fractions [17, 18|. This 
is in contrast with some different approaches proposed so far to accommodate 
the Coulomb force in three-nucleon (3N) continuum calculations: the screening 
and renormalization approach in momentum space Faddeev integral equations 



[19, 2C, 21, 22], the partial differential equation approach in configuration space 
[23, 24, 25], the three-potential formalism with the Coulomb-Sturmian separable 
expansion method f26|| , the Kohn variational method with the Pair-correlated 
Hyperspherical Harmonic (PHH) basis p7], 28]. 



Since the Coulomb-modified Faddeev equation was given heuristically in ref. 
fll2| , we present a formal derivation of the equation based on a time-independent 
scattering theory in Sect. |2[ which should be useful in further development in 
more general few-body problem. The modification is performed by introducing 
an auxiliary potential to cancel the long-rangeness of the Coulomb potential. 
As will be explained in Sect. ^, the cancellation works successfully at energies 
below the three-body breakup threshold. Numerical results at those energies are 
thereby given in Sect. |3| for realistic NN forces and models of 3NF. Summary is 
given in Sect. |j. 

2 Three-Body Scattering with Coulomb Force 

2.1 Notations 

We consider a system of three equivalent particles (nucleons) of the mass m. The 
Jacobi coordinates {x a ,y a } are defined as (see Fig. ||) 

x a = rp - r 7 , y a = r a - - (rp + r 7 ) , (2.1) 

where a, (3, and 7 denote 1, 2, and 3, respectively, or their cyclic permutations, 
and r a is the position vector of the particle a. The Hamiltonian of the system 
in the center-of-mass (cm.) frame is 

H = H + V = H + Vt + V 2 + V 3 , (2.2) 

where Hq is the internal kinetic energy operator of the three-body system, 

h 2 3h 2 

H = T x (x a )+T y (y a ) = --Vl a - — Vj a , (2.3) 

and V a is an interaction potential for the pair (/5, 7), which consists of a short- 
range nuclear potential v s a {x a ) and the Coulomb potential v^(x a ), 

V a = v s a (x a ) + v c a (x a ) = v s a (x a ) + Mj^, (2.4) 

Here, Zpe and Z^e are the electric charges of the particles (3 and 7, respectively. 
We are going to solve a 3N Schrodinger equation with the cm. energy E, 



HW = EW, 



(2.5) 
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Figure 1. The Jacobi coordinates of the three-body system. 



to obtain the scattering state & initiated by a state £\ , 

(x 1 ,y 1 \~ 1 )=<f> d (x 1 )F(y 1 ;p ), (2.6) 

where 4> d (xi) is the deuteron state function of the pair (2,3), and F(yi;po) is 
a free state function of the nucleon 1 with the momentum pq. These functions 
satisfy the following Schrodinger equations, 

[T x ( Xl ) + Vi(a!i)] = Ed<l> d (xi), (2.7) 

and 

3ft 2 

T y ( yi )F( yi ;p ) = E po F( yi ;p ) = —p 2 F( yi ;p ), (2.8) 

with 

E = E d + E P0 . (2.9) 
In a time-independent formal scattering theory, the scattering state is given 

as 

^ = lim ieG{E + ie)S 1 = lim &(e), (2.10) 

e^O e^O 

where G(z) is the total three-body Green's function, 
Using the free Green's function, Gq(z), 

z — -on 
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and the resolvent relation, 

G(z) = G (z) + G (z)VG(z), (2.13) 

we obtain 

9(e) = ieG (E + ie)S 1 + G (E + ie)V9(e). (2.14) 

In the Faddeev theory with short-range potentials, the three-body state is 
divided into three components (Faddeev components), 

9(e) = $ 1 (e)+$ 2 (e)+$ 3 (e), (2.15) 

where the Faddeev components are defined as 

9 a (e) = ieG (E + ie)E 1 5 aA + G (E + ie)V a 9(e). (2.16) 

Eqs. (|2~l5| ) and ( ^Wj j give 



$ a (e) = ieG a (E + ie)Si 6 aA + G a (E + ie)U a (^(e) + <P 7 (e)) , (2.17) 
where we have defined the channel Green's function G a (z), 

G a (z) ee J-—- = 1 G (z). (2-18) 

2 - ilo - Vq 1 - Go(z)Va 



By taking the e- limit in Eq. (|2.17 ), we reach the Faddeev equation 



<2> a = Jim^a(e) = Hi <W + G a (£ + iQ)K, (*/3 + * 7 ) • (2.19) 

However, the right hand side of Eq. ( p.!9| ) includes the following term, 

G (S + iO)<(^ + # 7 ), (2.20) 

which is known to cause a severe singularity due to the long-range Coulomb 
force, v^(x a ). 

2.2 The Coulomb-Modified Faddeev Equation 



The singularity in the kernel of the Faddeev equation (|2.19|) occurs due to the 
long-rangeness of the Coulomb interaction v^(x a ). We therefore consider to in- 
troduce an auxiliary potential that cancels the long-rangeness of v^(x a ) for a 
large x a . For this purpose, we define a Coulomb potential up (yp) that acts be- 
tween the cm. of the pair (7, a) and the spectator j3 with respect to the charges 
of the pair (/3, 7) (see Fig. ||), 

= (2-21) 

V0 



Using an expression, 



yp = x a + ~ X0 , (2.22) 
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which is easily given with the definition of the Jacobi coordinates, Eq. ( |2.1| ), we 
can show that 

111 1 0{x-% (2.23) 



C a \X a + 2*E/3\ 



when x a becomes large with keeping xp finite. The finiteness of xp is implicit 
in a Faddeev component <Pp, because the pair (7, a) is in the deuteron bound 
state or in a closed channel as far as we stay at energies below the three-body 
breakup threshold. (See the redefinition of the Faddeev component below, Eq. 
(ET27D.) The range of v^(x a ) — Up^(yp) is therefore short if it is combined with 
&p, i.e. 



\v c a (x a ) - u c p n (yp)}$p = ZpZ^e 2 I — - — J $p -» ZpZ^OixT^p. (2.24) 
1 ' \ x„ yp I 



d a 



y 




P 



Figure 2. Relation between the Coulomb potential Va(x a ) and the auxiliary Coulomb potential 
u% ~(yp), which acts between the spectator /3 and the cm. of the pair (7, a). 



With this idea, we write the second term on the right hand side of Eq. ( 2.14 ) 
as follows: 

G (E + iE)V&(e) = G (^ + ie)|^i*i(e) + ^2*2(e) + ^3(e)} 

+G (E + is) I (Vi - u c 2 , 3 ) <P 2 (e) + (Vx - u c 3j2 ) #3 (e) 



+G (E + k) j(V 2 - u c 3<1 ) 3 (e) + (V 2 - uf, 3 ) 
+G (£ + k)| (V 3 - uf i2 ) *i(e) + (V 3 - ^ 2 (e)} 
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+G (E + is)! [ u c 1 $ 1 {e) +y%$ 2 (e) + «§<P 3 (e)}, (2.25) 

where the last term on the right hand side including a potential u c a , which is 
defined as a two-body Coulomb interaction between the pair (/3, 7) and the spec- 
tator a, 

< = <(y*) = <AV«) + ««,>«) = Za ^ + Z ^ e \ (2.26) 

Vol 

is added to compensate the advisedly subtracted terms with ii^»'s. 

We correspondingly redefine the Faddeev components instead of Eq. (|2.16| ) 

as 

* a (e) = ieG (E + <5 a)1 + G (£ + le){V a + u£)<P a (e) 
+G (£ + ie) (V a - u% y ) ^(e) 

+G (£? + i£)(Va-i4 )i9 )* 7 (e), (2.27) 

which gives 

* a (e) = ie£ a (£ + te)£i <f Q ,i + + ie) (v a - u% y ) $ p (e) 

+g a (E + ie)(y a -ulp)$ 1 (e), (2.28) 

where we have defined the channel Coulomb Green's function Q a (z) as 

= w\i c = -1 r < uv I c , G o(4 ( 2 - 29 ) 

Taking the e-limit in Eq. ( |2.28[ ) and using a relation, 

limieixuy^iE + ie)~!) = ^(x^F^y^po,!]!), (2.30) 

e^O 

where F c (yi;pQ,r]i) is a scattering function in the Coulomb potential itf(yi), 
[^(yi)+«5(yi)]F c (y i; p ,?7i) = ^o^ c (l/i;Po^i), (2.31) 



and rj a is the Coulomb parameter, 



= S ^ + Z > 2 , (2.32) 
we get the Coulomb-modified Faddeev equation, 

$ a = <^F c (p , m) s a ,t + G a (E + iO) (y a - u c pn ) 

+G a (E + iO) (y a - u c ljP ) * 7 . (2.33) 

Since the potential term V^, — is a short-range interaction with respect 
to |a; a | as mentioned above, Eq. ( |2.33 ) can be treated as in the nd problem. At 



the same time, it is noted that the appearance of the channel Coulomb Green's 
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function Q a (E +i0) expresses that the spectator particle is distorted by the long- 
range Coulomb force. In the original idea of the Faddeev theory with a short- 
range interaction, while pair particles are interacting, the spectator particle is 
considered to be free. This picture is no more valid with the Coulomb force, 
and the auxiliary Coulomb potential creates the Coulomb wave function that 
expresses the distortion. 

For the case of identical particles, the initial state should be 

S = Si + S 2 + S 3 . (2.34) 

In this case, each Faddeev component has the same functional form and satisfies 
the same equation, 

<2> Q = ^(po^a) + Ga(E + iO) {[V a - ufrj&f, + [V a - u c ^]<P y } . (2.35) 
2.3 Elastic Singularity and Elastic Amplitude 

In this subsection, we derive an expression for the elastic amplitude for the pd 
scattering below three-body breakup threshold. For this, we write the Coulomb- 
modified Faddeev equation ( 2.35|) as 

# a = cf> d F c ( Po , m) + GiA^i, (2.36) 

where the operator Ai involves a coordinate transformation, 

= (y x - u c 2>3 ) <2> 2 + (Vi - u§ i2 ) # 3 - (2.37) 

Hereafter, we drop the particle number for simplicity. In the case of the 
scattering problem, the Green's function Q possesses a pole corresponding to the 
deuteron bound state. In order to treat this pole, we apply a subtraction method, 
in which we insert an identity, 

1 = |^)(^| + [l - |^)(^|] , (2.38) 

between Q and A in Eq. ( |2.36| ) to obtain 

(x,y\<P) = ( j) d (x)F c (y;p ,r ) ) + 4> d {x)(y\G c (4> d \A<I>) + (x,y\rA$). (2.39) 

Here, G c is the Coulomb Green's function for the outgoing particle, 

6c = -■ £TT> ( 2 ' 4 °) 

E po +ie-Ty(y) - u c {y) 



whose analytical form is known |29j as 

(y\G c \y') =J2yt m (y)G^(y,y')Yr(y' ), (2.41) 

and 

<&W) = -^ P0 e ^ ){P0y > )F ^ y <\ (2.42) 
3h 2 p y> p y< 
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where i*)(r) is the (real) regular Coulomb function, u^~\poy) is defined with the 
(real) irregular Coulomb function Ge(r) as 



,(+), 



*£ (poy) -- 

and a i is the Coulomb phase shift, 

a e = argr(£+ 1 + irj). 

The operator r is 



(2.43) 



(2.44) 



r = 



i 



E + ie - H - V(x 
dp\F c (p, V )) 



u c (y) 



l-\<f> d )(<t>° 



1 



\(/> d )(<p d \ 

E q + \e-T x {x)-V{x) ~ E g -E d 



where E„ is the energy of the two-body subsystem, 



E„ = E 



3^ 
4m 



P 



ti 2 



in 



(2.45) 
(2.46) 



and F c (p, rj) is the Coulomb-modified spectator state, which is a solution of Eq. 
( 2.31| ) with the eigenvalue of ^^p 2 - Because the singularities in the first and 
second terms on the right hand side in Eq. ( 2.45| ) cancel each other, we can 
apply the standard quadrature to perform the p-integration in Eq. (|2.45| ) if both 
terms are treated together. 

The first and second terms on the right hand side of Eq. ( 2.39| ) may be 
considered as a state that describes the elastic scattering. We correspondingly 
define a elastic function, J- (y), as 



T e \y) = F c (y; Po , v ) + (y\G c (<t> a \A$). 



(2.47) 



Substituting a formal solution of Eq. ( |2.39| ) for ^ to the second term on the right 
hand side of Eq. ( 2.47 ), we obtain the following equation, 



T e \y) = F c (y;p , v ) + (y\G c (<j> d \A 



i - rA 



b d T el ). 



(2.48) 



Although the interaction term in Eq. (2.48) is very complicated because of the 
coordinate transformation operator A, this equation has the same form as the 
Lippmann-Schwinger equation for a two-body scattering problem with a short- 
range force and the Coulomb force. As in the case of the two-body problem, 
using the explicit form of the Green's function G c , Eqs. (2.41) and (2.42), the 
asymptotic form of the function u^(j>oy), 

4 +) (Poy)\y-+oo = ex.p(+i(p y - r]ln2p y - e-Tr/2)) , (2.49) 
and the decomposition of the Coulomb scattering function F c (y;po,r]), 

,Fi{poy) 



F c (y-po,v)=^Y. iiY r*(P°) Y t m (y) e 

£,m 



poy 



(2.50) 
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we obtain an asymptotic form of the elastic wave function T e (y), 

T el (y)\ y ^ oo = F c (y;p ,7 1 ) + r l (y) , (2.51) 

y 

where f el (y) is defined as 

f e \y) = -^^(F-{p oy ,r,)^\A^). (2.52) 
The scattering amplitude of the elastic scattering is therefore given as 

f(y) = f c (y) + f el (y), (2.53) 

where f c (y) is the pure Coulomb scattering amplitude, 

fly) = * e -ivH S m*8/2)+2ia _ (2.54) 

y ' 2p o sin 2 0/2 v ; 

3 Numerical Results 

3.1 Some Remarks on Numerical Methods 

In this section, we present some results of our calculations to solve the Coulomb- 
modified Faddeev equation ( |2.36|) for the pd scattering problem at energies below 
the three-body breakup threshold. Although techniques used in the numerical 
calculations are essentially the same with those in the 3N bound state studies 
[14, 16, 17, pi], we give some remarks in this subsection. 



We solve the Coulomb-modified Faddeev equation fl2.36Q with an iterative 



method, which is called as the Method of Continued Fractions [[D], 18]. In most 
iterative methods to solve a linear integral equation, a A;-th order process involves 
the calculation of a function 17 given by operating the integral kernel to a known 
function Q that is produced from the (k — l)-th order solution. In the present 
case, we have 

n{x, y) = (x, y\QAQ) = <p d (x){y\G c (cP d \Af2) + (x, y\rAQ). (3.1) 
The first term on the right hand side of Eq. ( |3.1] ) is calculated numerically 



using the explicit expression of the Green's function G c , Eqs. (2.41) and ( |2.42j ). 

From the definition of the -T-operator, Eq. ( |2.45| ), the second term on the 
right hand side of Eq. ( |3.1[ ) reads 



(x,y\rA\f2) = J d P F c (y-p, V ){0(x;p)- ( l ) d (x)C(p)}, (3.2) 
where C(p) is given by 



C(p) = -(cf> d F c (p, V )\An), (3.3) 

&q — &d 



and &(x;p) by 
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In actual numerical calculations, this 0(x; p) is obtained as a solution of Schrodinger- 
like differential equation, 

[E q - T x {x) - V(x)] &(x;p) = (F c (p, rj)\Af2), (3.5) 

with the boundary condition that 0(x;p) vanishes at a large x because E q is 
negative in the present case. (See Eq. Q2.46p .) After a decomposition of 0(x;p) 
into partial waves, we have a set of ordinary differential equation, whose solutions 
can be obtained with the usual technique as in the two-body problem. 

For the partial-wave decomposition, we use the following angular-spin-isospin 
functions in the J — j coupling scheme, 



25+1 



Lj,£ f JSM ) 



where L is the relative orbital angular momentum of the pair nucleons; S the 
total spin of the pair; J the total angular momentum of the pair (J = L + S); 
i the orbital angular momentum of the spectator nucleon; j the sum of I and 
the spin of the spectator (J = £ + ~); Jo and Mq the total angular momentum 
of the three nucleons (Jo = J + j) an d its third component, respectively; tt 
the parity of the system (= (— ) L+e ); I and Mj the total isospin and its third 
component of the pair; /i the third component of the isospin of the spectator; 
Mt the third component of the total isospin of the three nucleons (Mt = +\ 
for the pd scattering and Mt = —\ for the nd scattering). Note that we use 
the isospin function in which the isotriplet pair nucleons are assigned as proton- 
proton (neutron-neutron) for Mi = +1(— 1) and proton-neutron for M; = 0. 
With this form, we take into account the mixing of the total isospin, 1/2 and 
3/2, arising from the inclusion of the Coulomb force and/or a charge symmetry 
breaking (CSB) NN forces. Actual form of the isospin functions, whose phases 
are chosen for convenience' sake in the coordinate transformation, is shown in 
Table | H Ell- 



in the present calculation, 3N partial-wave states for which NN force acts are 
restricted to those with total two-nucleon angular momenta J < 3 (see Subsect. 



3.2), and the total 3N angular momentum Jq is truncated at Jq = 19/2 in 



calculating scattering observables in Subsects. 13^ and 



Table 1. Three- nucleon isospin functions for Mt = ±i 



I Mj \i \IMj,^;M T ) 

#0 ±| ±|00)[|±i) 

#1 1 ±j -|io)|i±j) 

#2 1 ±1 T J |i±i)||t|) 
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3.2 Phase Shift Parameters 

As a first step, we use the Malfliet-Tjon (MT) I-III NN potential @ with pa- 
rameters given in ref . |3lf | . Because of the simplicity that the MT I-III potential 
acts only on the spin triplet and the spin singlet S-wave NN states, the 5"-wave 
nucleon-deuteron (Nd) elastic channel is classified only by the channel spin: the 
sum of the nucleon spin and the deuteron spin, namely, 1/2 (the doublet state) 
or 3/2 (the quartet state). In Table calculated S- wave phase shifts of the 
doublet ( 2 5) and the quartet ( 5) states for Nd elastic scattering at some en- 
ergies below the three-body breakup threshold are compared with the results 
of the configuration-space Faddeev calculations PH , whose results for a realis- 
tic NN potential model, the Argonne V14 model (AV14) |p2| , were published as 
benchmarks |33|] in the pd calculation. 



The differences between our calculations and those of ref. [23] are an order 
of 1% at every energies. 

Table 2. The nd and pd S-wave phase shifts with for the MT I-III potential. The phase shifts 
are in degrees. 



Elf (MeV) Ref. || This work Ref. g| This work 



0.0015 


-0.230 


-0.228 


-2.09 


-2.09 


0.15 


-3.28 


-3.26 


-20.4 


-20.5 


1.5 


-20.7 


-20.6 


-55.8 


-55.9 


2.45 


-28.6 


-28.6 


-66.7 


-66.7 


3.27 


-33.6 


-33.4 


-73.6 


-73.5 



E^ ao (MeV) Ref. |2J This work Ref. J23J This work 



0.15 


-0.537 


-0.527 


-7.46 


-7.44 


1.0 


-10.6 


-10.5 


-37.7 


-37.2 


1.5 


-16.2 


-16.2 


-46.5 


-46.4 


2.0 


-21.1 


-21.0 


-53.5 


-53.3 


3.0 


-28.8 


-28.7 


-63.8 


-63.4 



Next, we examine our calculations with a more realistic NN potential, the 
Argonne Vig model (AV1 8) [|34[ . We follow the convention of the phase shifts and 
mixing parameters in ref. |35[| , where the Nd elastic channel is denoted as 2K+1 £j 
with K being the channel spin (K = 1 + ^, Jo = K + £). These parameters 
are calculated with partial-wave amplitudes in the channel spin scheme, which 
are transformed from our partial-wave amplitudes in the J — j coupling scheme 
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in a trivial manner . In Tables [3| and ||, we present results of the phase shift 
parameters up to Jo = 5/2 at Ejf h = 1 MeV and E]f h = 3 MeV, respectively. In 
the tables, our calculations with J < 2(3) are denoted as J max = 2(3). The results 



are compared with those of the variational method with the PHH basis |27j, 



whose results for the AV14 potential were also published as the benchmarks 



[35, 33]. From these tables, we see that our calculations almost converge at 
Jmax = 3 and the differences from the PHH calculations are again order of 1%. 
In addition, although results are not shown in the tables, we have examined 
the convergence of the partial-wave expansion of the term v^,(x a ) — u% (yp) in 
Eq. ( |2.35| ) by neglecting J = 3 NN force components v s a {x a ) in the J max = 3 
calculation. The discrepancies between the results and the J max = 2 calculations 
are much less than 1%, which shows that the effect of the J = 3 components in 
the term v^(x a ) — u c p Ayp) is negligibly small. 

Tables ||, |3[ and [| demonstrate that our formalism gives quite promising 
results in calculating the nd and pd elastic scattering at low energies. 

Having established the accuracy of our calculations, we proceed to compare 
our results with available experimental analysis. In ref. 0], a phase shift analysis 
(PSA) was performed for the pd scattering at Ep ah = 3 MeV, corresponding to 
£Lab = 6 Me y_ Resu i ts f 

our calculations for the phase shift parameters are 
compared with the PSA results in Table [|. The authors in ref. Q noticed two 
significant differences of their PSA results from those of Faddeev calculations 



with the PEST16 NN potential 2C[: the 2 S , 1 / 2 phase shift and the S — D mixing 
parameter, f?i/2+- The similar differences are observed for the comparison with 
our AV18 results in Table [|. The discrepancy in the 2 5' 1 / 2 phase shift is related 
to the binding energy of the 3N bound state. While the experimental value of 
3 He binding energy is 7.72 MeV, the calculated value for AV18 is 6.91 MeV. 
To solve this discrepancy, we introduce a 3NF model by Brazil group (BR) 
based on the two-pion exchange process with the cut-off parameter of 700 MeV. 
Calculated 3 He binding energy for the AV18 plus the BR-3NF (AV18+BR) is 
7.79 MeV. At the same time, although not enough but a significant improvement 
is observed for the 2 S1/2 phase shift and the 771/2+ parameter as shown in Table 

In addition, we observe a significant discrepancy in the splittings of the spin 
quartet P-wave phase shifts, dip with Jo = 1/2, 3/2, and 5/2, although their 
average for AV18 and AV18+BR agrees with the PSA values. This difference is 
closely related with so called the A y {9) puzzle: calculated values of the nucleon 
analyzing power A y (9) with realistic NN potentials for energies below ~ 30 MeV 



underestimate systematically from experimental values [37, [38]]. In ref. 3J3], a 
combination of P-wave phase shifts, 

" ^p 1/2 ~ 5^P 3/2 + 9 ^P 5/2 = 4 (V 3/2 " ^ 1/2 ) + 9 (V 5/2 " ^P 3/2 ) » ( 3 - 7 ) 
was shown to be a good measure of the height of ^L-peak. The value of Eq. 



(P77D is 10.05±0.70 for the PSA values, 7.15 for AV18, and 7.88 for AV18+BR. 
The values of AV18 and AV18+BR are smaller than the experimental values by 
about 30%, which is consistent with the underestimation of the A y -peak height 
at 3 MeV. (See Subsect. O.) 
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Table 3. The nd and pd phase shifts parameters at E^f = 1 MeV with AV18 potential. The 
PHH results are taken from ref. |2q]. Values of all parameters are in degrees. 



nd pd 
This work PHH This work PHH 

Jq Jvoax — ^ </max — 3 ^max — 2 </max — 3 



i+ 4 n 

S\/2 
V1/2+ 

2 ^1/2 
4 A/2 
e l/2- 

3+ 4c , 
2 ^3/2 

-^3/2 

4 Av2 

£ 3/2+ 
C3/2+ 
%/2+ 

3- 4 77 
2 ^3/2 

^3/2 

4 ^3/2 

e 3/2- 

C3/2- 

%/2- 

5+ 4r< 

2 ^5/2 
2 A/2 
4 A/2 
£ 5/2+ 
C5/2+ 
%/2+ 

5- 4p 

F5/2 
4F 5/2 
e 5/2- 
C5/2- 
115/2- 



-0.99 
-18.5 
0.875 

-4.13 
12.1 
3.46 

-46.9 
0.572 
-1.06 
0.607 
0.516 
-0.110 

0.124 
-4.08 
14.0 
-1.24 
-0.173 
-1.05 

-0.015 
0.568 
-1.12 
-0.276 
-0.273 
-0.810 

13.3 
-0.064 
0.129 
0.441 
0.383 
-0.120 



-0.99 
-18.1 
0.919 

-4.15 
12.1 
3.44 

-46.9 
0.571 
-1.06 
0.606 
0.515 
-0.107 

0.124 
-4.09 
14.0 
-1.28 
-0.177 
-1.04 

-0.015 
0.567 
-1.12 
-0.265 
-0.270 
-0.818 

13.2 
-0.064 
0.129 
0.443 
0.390 
-0.123 



-0.98 
-18.1 
0.928 

-4.13 
12.0 
3.47 

-46.7 
0.564 
-1.05 
0.621 
0.511 
-0.107 

0.121 
-4.08 
13.9 
-1.24 
-0.177 
-1.04 

-0.015 

0.560 

-1.11 

-0.277 

-0.272 

-0.821 

13.2 
-0.063 
0.127 
0.447 
0.390 
-0.123 



-0.770 
-13.8 
0.93 

-3.32 
8.98 
2.96 

-36.9 
0.432 
-0.821 
0.628 
0.493 
-0.093 

0.095 
-3.29 
10.5 
-1.05 
-0.184 
-0.968 

-0.011 
0.429 
-0.869 
-0.272 
-0.272 
-0.818 

9.9 
-0.049 
0.099 
0.414 
0.385 
-0.120 



-0.770 
-13.5 
0.98 

-3.33 
8.95 
2.95 

-36.9 

0.431 

-0.823 

0.613 

0.493 

-0.091 

0.095 
-3.29 
10.4 
-1.08 
-0.187 
-0.961 

-0.012 
0.428 
-0.871 
-0.249 
-0.270 
-0.824 

9.8 
-0.049 
0.099 
0.419 
0.391 
-0.122 



-0.771 
-13.2 
1.02 

-3.34 
9.19 
2.89 

-37.0 
0.442 
-0.829 
0.784 
0.500 
-0.089 

0.097 
-3.30 
10.7 
-1.02 
-0.182 
-0.932 

-0.011 
0.439 
-0.879 
-0.364 
-0.300 
-0.943 

10.0 
-0.049 
0.100 
-0.218 
0.384 
-0.127 
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Table 4. The nd and pd phase shifts parameters at = 3 MeV with AV18 potential. The 

PHH results are taken from refs. [E^, E3. Values of all parameters are in degrees. 



nd pd 







This work 

<^max — 2 Jmax — 3 


PHH 


This work 

•Anax — 2 Jmax — 3 


PHH 


1 + 

n 
Z 


4 D 1/2 


-3.86 


-3.87 


-3.85 


-3.55 


-3.56 


-3.56 






-35.8 


-35.2 


-35.3 


-32.8 


-32.3 


-32.2 




771/2+ 


1.03 


1.10 


1.12 


1.02 


1.07 


1.10 


i - 

2 


1/2 


-7.42 


-7.47 


-7.49 


-7.32 


-7.37 


-7.36 




4 Pl/2 


24.3 


24.2 


24. 2 


21.8 


21.7 


22.1 




e l/2- 


6.70 


6.67 


6.68 


5.76 


5.74 


5.71 


3 + 
2 


4o 
^3/2 


-70.0 


-70.1 


-69.9 


-62.9 


-63.0 


-63.1 




2 n 

/ /•> /(-> 


2.39 


2.37 


2.36 


2.09 


2.07 


2.15 




4 -Dq/9 


-4.15 


-4.16 


-4.14 


-3.81 


-3.82 


-3.83 




^3/2+ 


0.724 


0.720 


0.747 


0.792 


0.753 


0.800 




C3/2+ 


1.37 


1.36 


1.35 


1.29 


1.29 


1.30 




%/2+ 


-0.377 


-0.363 


-0.363 


-0.327 


-0.320 


-0.322 


3 - 
2 


4 7? 


0.931 


0.929 


0.920 


0.830 


0.828 


0.849 




- 1 3/2 


-7.10 


-7.12 


-7.18 


-7.08 


-7.11 


-7.14 




^ Of 2 


26.1 


25.9 


26.0 


24.0 


23.9 


24.2 




°d/z — 


-2.63 


-2.69 


-2.61 


-2.24 


-2.30 


-2.20 




C3/2- 


-0.245 


-0.261 


-0.265 


-0.305 


-0.317 


-0.321 




%/2- 


-3.57 


-3.52 


-3.52 


-3.14 


-3.10 


-3.11 


5 + 
2 


4/^f 
<-*5/2 


-0.208 


-0.208 


-0.206 


-0.189 


-0.189 


-0.189 




2 A/2 


2.36 


2.34 


2.33 


2.07 


2.05 


2.13 




4 A/2 


-4.48 


-4.49 


-4.46 


-4.09 


-4.11 


-4.13 




£ 5/2+ 


-0.316 


-0.285 


-0.315 


-0.315 


-0.257 


-0.350 




C5/2+ 


-0.716 


-0.693 


-0.701 


-0.707 


-0.688 


-0.699 




%/2+ 


-1.97 


-2.03 


-2.04 


-2.04 


-2.09 


-2.07 


5- 
2 


^5/2 


26.2 


26.1 


26.0 


23.8 


23.7 


23.9 




^5/2 


-0.467 


-0.470 


-0.466 


-0.430 


-0.434 


-0.433 




4 -^5/2 


0.961 


0.957 


0.951 


0.862 


0.858 


0.876 




e 5/2- 


0.481 


0.495 


0.538 


0.412 


0.434 


0.343 




C5/2- 


0.898 


0.927 


0.926 


0.897 


0.928 


0.932 




%/2- 


-0.325 


-0.340 


-0.334 


-0.338 


-0.349 


-0.343 
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Table 5. Phase shift parameters for pd scattering at E p a =3 MeV calculated for AV18, 
AV18+BR and AV18+BR+LS1. Values of PSA are taken from ref. Q. Values of all parameters 
are in degrees. 



PSA 



AV18 AV18+BR AV18+BR+LS1 



-63.95±0.28 


-63.01 


-63.02 


-63.02 


OA 87-1-0 35 


39 9Q 

OZ . & u 


97 3^ 
z / .oo 


97 43 

& I .4:0 


23.37±0.11 


23.09 


23.29 


23.19 


n m _i_n nfi 


n 1 7 

u. 1 ( 


n D4 


n 4n 


9 AQ-l-n 1 1 


9 17 


1 88 
1 .oo 


9 1 8 

Z. lO 


7 1 1 _i_n OA 


7 94 


7 94 


7 9fi 


-0.14±0.25 


0.26 


0.29 


0.34 


o. / ^imu.ui 


O.Ol 


3 81 

O.Ol 


3 81 

O.Ol 


n qi _i_n hq 
U.ol±U.Uo 


n Q7 
U.o 


n Q7 


U.oo 


-0.28±0.03 


-0.29 


-0.29 


-0.28 


-0.16±0.03 


-0.26 


-0.27 


-0.27 


1.99±0.09 


2.06 


2.06 


2.06 


0.0 (fixed) 


-0.02 


-0.02 


-0.02 


2.00±0.10 


1.07 


1.51 


1.50 


1.20±0.04 


1.29 


1.29 


1.29 


-0.28±0.03 


-0.32 


-0.31 


-0.31 


5.73±0.13 


5.74 


5.78 


5.69 


-2.47±0.05 


-2.30 


-2.41 


-2.52 


-1.25±0.54 


-0.32 


-0.31 


-0.30 


-3.13±0.23 


-3.10 


-3.11 


-3.12 


-0.62±0.20 


0.93 


0.92 


0.91 


-0.45±0.04 


-0.35 


-0.35 


-0.35 


2.15±0.24 


0.75 


0.76 


0.78 


-0.75±0.12 


-0.26 


-0.26 


-0.25 



S-wave phases: 

£3/2 

2 S 1/2 
P-wave phases: 

4 Pj average 
P5/2 — A P%/2 
P3/2 — A P\/2 

2 Pj average 

2 P3/2 — 2 P\/2 
I?- wave phases: 
A Dj average 
A D 7/2 - 4 D 5/2 

4 D 5/2 -*D : 
^3/2 



4 -Dq/9 — A P > 1/ 



3/2 
2 



2 Dj average 

2 £5/2 - 2 Av2 

S — D mixing parameters: 
V1/2+ 
C3/2+ 
V3/2+ 

P — P mixing parameters: 
£l/2- 

£ 3/2- 

P — F mixing parameters: 

C3/2- 
V3/2- 
C5/2- 
V5/2- 

D — D mixing parameters: 

£ 3/2+ 
£ 5/2+ 
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Although we do not have any consensus to solve the A y (0) puzzle at the mo- 
ment, one possible reason for the discrepancy might be an unknown three-nucleon 



force that has a character of spin-orbit forces. In ref. [40], Kievsky introduced 
purely phenomenological spin-orbit three-nucleon potentials. We take one of the 
3NF models that has the longest range (LSI). The binding energy of the 3 He for 
the LS1-3NF in addition to the BR-3NF (AV18+BR+LS1) is 7.77 MeV. Results 
of the phase shift parameters for AV18+BR-I-LS1 are also shown in Table || The 
LS1-3NF gives a remarkable effect on the phase shift splittings of 4 -Pj states. 



As a result, the value of Eq. (gj) for AV18+BR+LS1 is 12.32. 



3.3 Differential Cross Sections 

In ref. ||, the angular distributions of the differential cross section (DCS) for the 
pd elastic scattering were measured precisely in the energy range of 2 < Z?p ab < 
18 MeV. In Fig. |, we plot their results at E^ ah = 3.0 MeV being normalized to 
our theoretical results of the AV18 calculation. We observe that the theoretical 
results overestimate the experimental data by about 5% at all scattering angles. 
Similar discrepancy has been observed in recent pd experimental data at low 



energies M 11]. 




Figure 3. Calculated and experimental pd differential cross sections normalized by the result 
of the AV18 calculation at i£ ab = 3.0 MeV. Dotted curve denote the result of the approximated 
Coulomb treatment with the AV18 nd calculation. Solid and dashed curves are results with the 
AV14', and dTRS' calculations, respectively. Experimental data are taken from ref. |jj]. 

The authors in ref. || extracted the cross-section minimum, which appears 
around 6 c .m. = 100°, from their data, and compared their results with the predic- 
tion of Faddeev calculations with the Paris NN potential, in which the pd scat- 
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tering amplitudes are approximated as the sum of the pd two-body Rutherford 
amplitudes and nd Faddeev amplitudes with the pd two-body Coulomb phases 
as was done in ref. [41]. As a result, they found that the ratio of the theoretical 
cross-section minimum to the experimental one reveals a strong energy depen- 
dence. The theoretical cross-section minimum overestimates the experimental 
one by about 20% at Ep ah = 2 MeV. The overestimation reduces to about 10% 
at 3 MeV, and once disappears at 5 MeV. At higher energies, the calculations 
underestimate the experimental one. (See Fig. 17 in ref. |||.) The dotted line 
in Fig. H|, which denotes the result of AV18 with the same approximation as 
that in ref. [41], is consistent with the result of ref. ||. We therefore conclude 
that the overestimation of DCS at low energies observed in ref. [|| is due to the 
improper treatment of the Coulomb force, and it is reduced to about 5% by a 
proper treatment of the Coulomb force. 

Besides the AV18 model, we examine some other realistic NN interaction 
models: a super soft-core model of de Tourreil et al. (dTRS) jE| and the earlier 
version of the Argonne model, AV14 |3^]. Since the isospin-triplet component 
of these potentials was made as either proton-proton or neutron-proton force, a 
charge dependence is introduced so as to explain the difference among proton- 
proton, neutron-proton, and neutron- neutron 1 Sq scattering lengths [16]. The 
results with these models (AV14' and dTRS') are plotted in Fig. ||[ The figure 
shows that the differences of the DCS due to the different input NN potentials 
observed at c . m . ~ 60° and <9 c . m . ~ 180° are at most 3%. At <9 c . m . ~ 120°, 
where the dependence on the input NN potentials is small, the calculated results 
overshoot the experimental data by about 5%. 

As is seen in the results of the pd S- wave phase shift parameters, we may 
expect that the introduction of the BR-3NF improves the DCS. The results of 
calculations with the AV18+BR normalized to those with the AV18 for Ep 3 ^ = 
0.4 — 3.0 MeV are shown as solid curves in Fig. |] together with experimental 
results of refs. Qg. The fi gures show that the introduction of the BR-3NF im- 
proves remarkably the fit to the empirical cross sections. It is well known that the 
3N binding energy and the doublet scattering length are strongly correlated. The 
improvement of the differential cross section due to the BR-3NF therefore might 
be nothing more than that to reproduce the 3N binding energy. To demonstrate 
this idea, we examine a fictitious spin-independent 3NF of the Gaussian form 
(GS-3NF), 

V G = V G e W {-(—) 2 ~ (— ) 2 } + (c.p.). (3.8) 
r G r G 

Values of the parameters, which are determined so as to reproduce the empirical 
3N binding energy, are tq = 1.0 fm and Vq = — 40 MeV. The results with the 
GS-3NF (AV18+GS) are shown as dashed curves in Fig. |. It is noted that the 
purely phenomenological 3NF achieves the similar improvement. 



3.4 Nucleoli Vector Analyzing Powers 

In Fig. H], we show results of the proton analyzing power A y (9) at several ener- 
gies below three-body breakup threshold comparing with available experimental 
data HI- The results with the AV18 potential (dashed curves) undershoot the 



18 



Low Energy Proton-Deuteron Scattering 



a(6) / a (e) 

v ' AV18 V ' 



E Lab = 0.4 MeV 

p 


E Lab = 0.9 MeV fc 

' Jm k 

I* . — 




\ 

— -"si 


E Lab = 0.5 MeV 

p 


E p Lab =1.0MeV j& 


\ ••••••• 


• 

I 


E Lab = 0.6 MeV 

"... 

— • 

— 


E p Lab = 2.0 MeV 








E Lab = 0.7 MeV 

* 


E Lab = 2.5 MeV 

p 


\ •••• 

• • 




E Lab = 0.8 MeV tM 


E Lab = 3.0 MeV 

p 









30 60 90 120 150 30 60 90 120 150 

8 (deg) e 



Figure 4. The pd differential cross section normalized to those calculated with AV18. Solid 
curves are those with AV18+BR. Dashed curves are those with AV18+GS. Experimental data 
are taken from ref. i] for E^ b = 0.4 - 1.0 MeV, and from ref. j| for ££ ab = 2.0 - 3.0 MeV. 



experimental values of A y {9) (The A y puzzle). The inclusion of the BR-3NF into 
the nuclear Hamiltonian increases only a small amount of A y {6) as denoted by 
the solid curves. On the other hand, the phenomenological spin-orbit 3NF, LSI 
[ 40 1 , improves A y {9) considerably as denoted by dotted curves. It is noted that 
the small amount of overestimation in A y {9) with AV18+BR+LS1 above 2 MeV 
is because of the combination of the LS1-3NF with the BR-3NF since the spin- 
orbit 3NFs in ref. |4C] are made so as to reproduce the experimental values of 
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Ay (9) by themselves, 

A merit to accommodate the Coulomb interaction in 3N system is that we 
can discuss the effect of CSB with high quality by comparing calculations with 
and without the Coulomb interaction. In ref. [16|, CSB effects are studied for 3N 
bound states. It was demonstrated that the calculated values of the S H binding 
energy (B( 3 H)) and those of s He (Bi^He)) for various combinations of NN 
forces and 3NFs have a linear correlation. When we plot the results of B( 3 He) 
without any CSB effect (but only with the Coulomb force) against those of 
B( 3 H), these points are fitted by a line, which does not pass through the point of 
the experimental values. When we include a CSB NN force, which reproduces the 
experimental difference of the proton-proton and neutron-neutron 1 So scattering 
lengths, Aa = 1.5 fm, together with the higher order electromagnetic interaction, 
we can draw a line that shifts to pass the experimental point. The CSB force 
used in the study of the 3N bound states is a central one. On the other hand, 
information of non-central spin-dependent CSB forces is expected to be given 
from the study of polarization observables in the nd and the pd scatterings. So 
far, the nucleon vector analyzing power A y (9) was measured for the nd and the 
pd elastic scattering for some energies [f|, |44}| . The measured values of A y {9) 
have peaks at 9 C , m ~ 105° at low energies. The difference between the nd and pd 
Ay (9) peak heights, which depends weakly on incident energies below Ejf h = 10 
MeV, and is roughly 0.01. 




Figure 5. Analyzing powers of the incident proton in the pd scattering at _E p a = 0.8 — 3.0 
MeV. Dashed curves are calculations with AV18, solid curves those with AV18+BR, and dotted 
curves those with AV18+BR+LS1. Experimental date are from ref. Ill for Ep ab = 0.8, 0.9, and 
1.0 MeV, and from ref. § for E^ h = 2.0, 2.5, and 3.0 MeV. 
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In Fig. ||, the calculated values of the difference, A™ d — AP d \g c m =105° at E]f h 
= 3 MeV are plotted against their average values, \Ay d + A^ d ^j /2\g c m =105° • In 
this figure, the results with AV14' and AVT4'+BR+LS1 are plotted as open 
circles, those with dTRS' and dTRS'+BR+LSl as open triangles, and those 
with AV18, AV18+BR, and AV18+BR+3NF as open squares. Apparently, the 
results are divided into two groups that are linearly correlated: one is the results 
with AV14' or dTRS' NN potentials as guided by the dotted line, and the other 
is those with AV18 as guided by the dashed line. 

We remark that the difference between these lines is attributed to higher order 
electromagnetic (EM) interactions in the AV18 potential. The EM interactions 
include a spin-orbit force arising from the magnetic moment interaction, which 
contributes as a CSB force [34]. To demonstrate the effect of the EM spin-orbit 
CSB interaction (EMls), a result of calculation with AVI8 neglecting EMls is 
plotted as a filled square, and that with AVT4' including EMls as a filled circle 
in Fig. |6|. The former point almost lies on the dotted line, and the latter on 
the dashed one. We therefore conclude that the higher order electromagnetic 
interaction in the AV18 model explains the difference of the two lines in Fig. |(| 
and gives about 10% effect in the nd-pd analyzing power difference. 

In ref. [45], the effect of the EM spin-orbit force is reported to be small in Nd 
analyzing power. However, the deuteron was treated as a single spin-1 particle 
in ref. KfJ. The present result therefore shows an importance of a three-body 
treatment of the pd system. 
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Figure 6. A™ d - A vd vs. (A" d + A v y d )/2 at 6» c . m . = 105°. Experimental date are taken from 
kl, h3J. The meaning of the points is described in the text. 
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4 Summary 

We have formally derived a modified Faddeev equation to accommodate the long- 
range Coulomb interaction in three-nucleon (pd) scattering, which was originally 
presented in ref. |TI|. In the modification, we introduced an auxiliary Coulomb 
force that acts between a spectator proton and the cm. of the rest pair nucleons 
(proton and neutron). This interaction takes into account two issues: (1) the 
Coulomb distortion of the spectator particle and (2) the cancellation of the long- 
rang part of the original Coulomb force in the integral kernel for energies below 
the three-body breakup threshold. The modified Faddeev equation can be treated 
similarly as in the case of the nd scattering. 

We presented some numerical results to solve the Coulomb modified Faddeev 
equation for the pd elastic scattering at energies below the three-body breakup 
threshold for some realistic NN potential models without or with three-nucleon 
forces. Results of the phase shift parameters agree with the available accurate 
results by different methods. 

Differential cross sections for an NN potential overestimate the experimental 
data by about 5% at low energies. This overestimation is recovered by introducing 
three-nucleon forces that reproduce the 3N binding energy. 

CSB effects are studied for the nucleon analyzing powers in pd scattering 
and nd scattering at 3 MeV. We found that a spin-orbit force due to the nucleon 
magnetic moment interaction plays a significant role in the difference. 
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